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Abstract We study optical activity induced by curvature. The optical activity 
model we present has two phenomenological gyration parameters, within which we analyze 
three model cases, namely, an exactly integrable model, the Landau-Lifshitz model and 
the Fedorov model, these latter two are solved in the short wavelength approximation. 
The model background is a gravitational pp-wave. The solutions show that the optical 
activity induced by curvature leads to Faraday rotation. 



1 



1 Introduction 



The simplest instances of the phenomenon of optical activity, or gyrotropy, are 
(i) the Faraday rotation, i.e., the rotation of the polarization plane for a linearly 
polarized electromagnetic wave propagating in a medium, and (ii) birefringence, 
for left-hand and right-hand circularly polarized electromagnetic waves [1]. Optical 
activity appears in materials which show linear spatial dispersion of light, as in 
the case of some crystals [1]. Indeed, of the 32 symmetry classes of crystals only 
15 admit optical activity, and have been studied theoretically and experimentally 
[2]. Phenomenological models of optical activity [3]-[9] are of great importance in 
applications to optics [10]. 

Following Landau and Lifshitz [11], the gravitational field in vacuum can be 
considered as a specific sort of medium. An example of this, and as first shown by 
Drummond and Hathrell [12], vacuum interacting with curvature produces birefrin- 
gence for an electromagnetic wave propagating in a curved spacetime. A glamourous 
consequence, is the gravitational rainbow considered by Lafrance and Myers [13]. 
Other effects related to quantum electrodynamics in background fields (either grav- 
itational or electromagnetic) show that the propagation of light depends on various 
factors, such as, photon polarization, photon frequency (temporal dispersion), direc- 
tion of propagation, and spacetime location [14]-[16]. As has been accepted recently 
(see, e.g. [17] and references therein), the one-loop gravitationally corrected Maxwell 
equations appear to be mathematically equivalent to the classical electro dynamic 
equations for the model of a linear, anisotropic, nonhomogeneous media. In this 
context the influence of the curvature on the electromagnetic waves has been de- 
scribed in terms of spacetime dependent electric and magnetic permeability tensors 
(see [17] and [18]-[20]). 

In this paper, we focus on the phenomenon of optical activity induced by the 
curvature of spacetime. First, our description of optical activity is phenomenological. 
The interaction constants, i.e., the gyration parameters, considered in our model are 
not derived from quantum electrodynamics. Nevertheless, the approach is motivated 
by the same principles, related to quantum corrections to Maxwell equations in a 
curved background [14]- [16]. The modeling of optical activity is of interest not only 
from a microscopic point of view, but also as a phenomenological set up [21]. Second, 
we choose a pp-wave as the background spacetime. It is a background well suited 
for our purposes, since it is an exact solution of Einstein field equations and has an 
explicit time dependence. It has attracted much attention in the recent literature 
in connection to the fact that it is the Penrose limit of any spacetime [22]. The 
main goal of the paper is to present an exact solvable model describing the optical 
activity induced by curvature, along the line of papers [18]- [20]. The results show 
that the interaction of the electromagnetic field with the background curvature (a 
pp-wave here) produces a Faraday rotation of the polarization vector, as in the 
classical theory of optical activity. The rotation frequency happens to be a function 
of the electromagnetic wave frequency and the Riemann tensor. 
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2 Two-parameter model of optical activity 



2.1 The three-dimensional classical prototype 

Maxwell equations in a dielectric medium with vanishing current of free charges are 

f)R rtf) 
V-D = 0, V-B = 0, Vx£=-— . Vx/?=^. (1) 

These equations include as usual the three-vectors D and B (the electric and mag- 
netic induction vectors, respectively), as well as the three- vectors E and H (the 
electric and magnetic fields, respectively). V is the standard three-dimensional dif- 
ferential vector operator. In the framework of linear electrodynamics, Landau and 

— * — * 

Lifshitz [1] connect the D and E three-vectors as follows, 

D" = e'^'^Eb + Y'^'dbE^ , (2) 

where e"^ is the electric permeability tensor, 7"^'^ is the optical activity tensor of rank 
three (it describes the gyration properties of the material), and dh denotes partial 
derivative. In this subsection latin indices a, b, c arc Euclidean indices running from 
1 to 3. When the medium is isotropic or has cubic symmetry, the 7"^^ tensor reduces 
to a pseudo-scalar / 

UJ 

where a; is a constant with units of frequency, and e"*'^ is the three-dimensional 
totally skew Levi-Civita symbol [1] . The pseudo-scalar / is the gyration parameter 
(it has no units). Thus, for the case of isotropic medium the relation (2) has the 
three- vector form 

D^eE + —V X E, (4) 

U) 

where the standard definition of the vector product V x E = e dbEc has been 
used, as well as the representation of the electric permeability tensor Sab — ^^ahi 
which is valid for the isotropic case. 5ab is a Kronecker delta with Euclidean indices. 

When one puts the standard condition — fi°'''Hh, where fi"-^ is the magnetic 
permeability, one obtains the so-called non-symmetric Landau-Lifshitz's model for 
optical activity (see, e.g. [9]). More generally, one can represent the dependence of 
B in H by the equation, (see, e.g. [5, 9]) 

B-^fi-'Hb + r''dbH,, (5) 

where the tensor 7"*"^ describes the magnetic aspects of the optical activity. Analo- 
gously to the relation (4), in an isotropic medium one obtains 

B^fxH+^VxH, (6) 

OJ 
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where /i is the magnetic permeabihty of the isotropic medium and / is a new pa- 
rameter, playing the same role for the magnetic field as / plays for the electric one. 
Finally, using the simplest model for the description of the properties of optically 
active media (see, e.g. [5, 10]), one can replace e and by their values in vacuum 
{s — 1 and /I = 1). Introducing two new gyration parameters a and b, defined by 
a = — and 6 = — ^, we obtain the three-dimensional form of the two-parameter 
model describing the optical activity, 

3 = E + aV X E, B ^ H -hV X H . (7) 

It will be better, for what follows, when we discuss the covariant generalization of 
equations (7), to have H as a function of B. Thus, inverting the second equation in 
(7) and keeping only up to first order terms in b one obtains 

D ^ E + aV X E, H ^ B + bV X B . (8) 

This two-parameter set of equations has three important one-parameter special 
cases. For a — b one obtains a model showing optical activity, which is fully intc- 
grable in the covariant analysis in a gravitational pp-wave background, (see section 
3.3.1). For 6 = and a = —b one obtains models not fully integrable but which 
can be dealt in the short wavelength approximation. The 6 = model is called the 
Landau- Lifshitz model (see equation (2) and paragraph before equation (5), see also 
the first part of section 3.3.2). The a = —b model, or Fedorov model, is very spe- 
cial indeed, since the constitutive relations preserve the duality invariance between 
{E, H) and {D, B) already present in the Maxwell equations in the non-conducting 
and non-charged media when one also includes magnetic monopole source terms [23] 
(see the second part of section 3.3.2). 



2.2 Covariant representation of the three-dimensional clas- 
sical prototype in Minkowski spacetime 

If we want to have a covariant model of optical activity we should rely on the 
covariant equations of electrodynamics of continuous media [1], 

Vfc/f*'^ = 0, (9) 

VfcF*^*^ = , (10) 

where Vfc is the covariant derivative, H^^ is the electric-magnetic induction tensor, 
and F***^ is the tensor dual to the Maxwell tensor Emm i-e.. 



?*ik Akls 



2v^ 



Fis. (11) 



Here, y=e**^''' is the four-dimensional Levi-Civita tensor and e*'^'* is the completely 
anti-symmetric four-dimensional symbol with e°^^^ = —£0123 = 1- Latin indices 
i, J, ... run from to 3. 
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The constitutive equations, connecting the induction tensor W'^ with the Max- 
well tensor F^nn can be written in the following form 

jjik ^ C'^ranp^^^ ^ V, , (12) 

where the so-called material tensor C'*'^™" describes the properties of the linear 
response (it includes the electric and magnetic permeabilities, as well as the magneto- 
electric coefficients), and the tensor includes the coefficients of the linear 
dispersion. These tensors have similar symmetry properties with respect to indices 
transposition, namely, 

^ikmn ^kimn ^iknm ^mnik (13) 

j-^iksmn j-^kismn j-^iksnm (14) 

Nevertheless, the D**^*"*" tensor is not in general symmetric with respect to the 
transposition of the pairs (ik) and (mn). The last term on the right hand side 
of equation (12) is a covariant generalization of the three-dimensional constitutive 
equations (2) and (5). 

In order to construct exphcitly the tensor D*'^*"*" for the two-parameter model 
of the optical activity we recall, first, that the electric induction, magnetic field, 
electric field and magnetic induction four- vectors are covariantly defined as 

D' = H'^Uk , H' = H*'^Uk , = F'^'Uk , B' = F*'^Uk , (15) 

respectively. These vectors are spacelike and orthogonal to the velocity four- vector 

D% = = E%, H% = = B% . (16) 

For W timelike and normalized to unity, one can represent the Maxwell tensor in 
the form 

jItts 



Fmn = ^mn EpUg — ejnnlsV—9 B ^ E^Un — EJJ^ — Tj^nlB , (17) 

where is the generalized 4-indices Kronecker delta, and rj^ni is an anti-symmetric 
tensor orthogonal to C/' defined as 

Vmnl = V^emnlsU\ r^^'^i = e'^'' U s . (18) 

The tensor rjmni is a covariant generalization of three dimensional Levi-Civita symbol 
gabc rpj^g identity 

\'''%lra = -5lsUjU^ (19) 



2 

will be important in what follows 
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Analogously to the decomposition (17) one can represent the induction 
tensor ifj^ in the form 

= C DpU, - ems^ H'W = DiUk - D^U, - rj^H' . (20) 

Taking into account the following covariant form of the relations (8) 

D'^E' + a rf^^dmEi , W^B' + b rf^^dmBi , (21) 

as well as using the definitions (15), we obtain from (17) and (20) the covariant 
constitutive equations for the two-parameter model 

^ik ^ pik ^ ) ^^^p^s ^ ^^p.] I _ (22) 

Comparing (12) and (22) one sees that the D^'^'^"^"' tensor takes the explicit form 

1 f fqjmn 



2 [V=9 
Here A^* = g^^ — U^U^ is a projector. 

3 Curvature generalization of the two-parameter 
model of optical activity 

3.1 Generalization of the model 

The curvature generalization of the model of optical activity can be obtained by 
replacing in the relation (22) the Levi-Civita tensor e*-''"'*/ by the dual Riemann 
tensor 

1 Ismn 

j^*qjmn = _j^qj^^ __ _ (24) 

This replacement can be motivated as follows. Hehl and Obukhov [24] state that 
a good linear ansatz for the constitutive relation between the induction tensor Hij 
and the Maxwell tensor Fki is 

mnkl 



Hij — Fij + ^ijmnX"^"' , (25) 

where ■^^'^^^ is a constitutive tensor with the same symmetries as the Riemann ten- 
sor. Now, the work of Drummond and Hathrell [12], on the corrections to quantum 
electrodynamics induced by curvature, tells us that a natural candidate to ■)^"^"-^^ is 
precisely the Riemann tensor BJ^'^^^, Thus we can replace one for the other in (25) 
and obtain 

Hij^Fij + xB*j'''Fki, (26) 

where x is some phenomenological coefficient. Let us summarize the following ideas 
discussed above: (i) the idea that the one-loop gravitationally corrected Maxwell 
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equations are mathematically equivalent to the classical electrodynamic equations 
in the presence of a linear, anisotropic, non-homogeneous media, is a generally ac- 
cepted idea; (ii) the phenomenological approach to the modifications of Maxwell 
equations plays an important auxiliary role in an analysis of the efi^ects induced by 
curvature; (iii) when the spacetime curvature influences the photons as a sort of 
anisotropic, non-homogeneous, non-stationary medium, one can expect that there 
exist the analogues of the corresponding phenomena in classical electrodynamics, 
indeed, the optical activity, as a well-known effect of linear dispersion, has to be 
considered among them; (iv) formally, the models describing the effects induced 
by curvature can be formulated by the following way: the tensor coefficients from 
the classical constitutive equations have to be replaced by the tensor coefficients 
containing the Riemann tensor. 

Following this line, we can take the covariant constitutive equations (22), de- 
scribing the classical optical activity, and replace e^j™"^/ by the dual Riemann 
tensor Rnjmn ^pj^jg procedure yields 

jjik ^ jpik ^ 1 ^ifc ^j^*qjmn \^aUW + h/^^'] } , (27) 

where now the gyration parameters a and h have been replaced by new phenomeno- 
logical curvature gyration parameters a and 6, respectively The units of a and h are 
a ■ length^ and h ■ length^. To simplify the notation we drop the bars from now on 
on a and h. Now the Maxwell equation (9) takes the form 

^^p^k ^ \^^k^ y ^ j^^^i-n ^ ,p^^ y^^vijs + ^^p.] | ^ Q . (28) 

As usually, we represent the Maxwell tensor Fjfe by 

Fiu = Vi Afe - VfeA , (29) 

where Ai is the four- vector potential. With this definition equation (10) is trivially 
satisfied. 



3.2 Equations for optical activity induced by a gravitational 
pp-wave 

The previous sections contain the description of the model for an arbitrary spacetime 
metric. It is of interest to study this model in a particular background. The pp-wave 
solution is a background well suited for our purposes, and incidentahy has attracted 
much attention in the literature, notably in the gravitational fields generated within 
string theory, due to the fact that it is the Penrose limit of any spacetime [22]. The 
pp-wave, a gravitational wave, is an exact solution of General Relativity and has an 
explicit time dependence. The pp-wave is described by the metric [25] 

ds^ = 2dudv - \e''>'{dxY + e-^^{dx^y] , (30) 
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where 



ct — 



ct + 



u — 



V — 



are the retarded and the advanced time, respectively. The functions L and (3 depend 
on u only, and are coupled by the equation 



L" + L iP' 



0. 



(32) 



A prime denotes generically the derivative of the function with respect to its own 
argument, in particular here the derivatives of L{u) and with respect to the 
retarded time u. One can assume as an arbitrary function of u and then solve 
for L. The curvature tensor has two non-zero components 



Rt 



L^0 



(33) 



'■u2u ^■u'iu 

As for the dual Riemann tensor, it has also two nonvanishing components 

-^2u3« = ^Zu2u = ^\^u ■ (34) 

The Ricci tensor Ri^. and the curvature scalar R are equal to zero. 

For a medium (or an observer) at rest in the chosen frame of reference one has 

1 



= {s: + SI) 



V2 



Calculating the components of the induction tensor Hik (27) one obtains 

H = F 

n.l) -*- 11.1} 1 



(35) 



(36) 



P I a+b f) ( flP7 F 1 _L 



+ 



a—b 
2V2 



8 ( R* n^^ F \ R* n^"' F — R* nP^ F^ 



(37) 



I a—b 



d (r* af^F ^ -^-T^ R* 0^'^ F —T^B* 0^^ F- 



23 



-^23 + 



2^/2" 



-rip 



R\u-fJ^2Fpu — R^.u-iJ^I'Fl 



2«7M^ 



pv 



(38) 
(39) 



where Greek indices take the values 2 and 3 only, and Y\i arc the Christoffel symbols. 
These Hik satisfy the following equations (sec (9) for the metric (30)) 



dvHyu + g^'^d^Hy^ = , 



(40) 
(41) 
(42) 



Let us focus now on specific exact solutions of the Maxwell equations given by 
(36)-(42). 
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3.3 Solutions to the modified Maxwell equations 

Consider electromagnetic waves traveling along the axis. For such a case the 
general solution of the electromagnetic potential Ai can be written as 



A2 = A2{u,v), A3 = A3{u,v), A^ = 0, A^ = . 
This means that the non-zero components of the Maxwell tensor are, 

■^Ua 9^Aq, , OyA^ . 



(43) 



(44) 



For this case the equations (40), (41) become trivial. Introducing the new functions 
M2{u,v) and M2,{u,v) by 



^2 = e^Ma (m, v) , A2 = e'^M^ («, v) , 



(45) 



and after some manipulations, one can reduce the two components in (42) to the 
two following coupled equations 



dud, [M2 + fjl 9„ {Rl^uM^) 



+ 



a~b 
4V2 



{Rl^uM,) =0 , (46) 



M3 + f±| a„ {R\,uM2) 



+ 



a—b 
4V2 



{Rl,uM2) =0 . (47) 



Note that the equation (47) can be obtained from (46) by interchanging the indices 
2 and 3, as well as 13 with —j3. 

Now, a simple model is obtained when a = b. In this case, the second part in 
equations (46) and (47) are zero, one has a D'Alembertian operator alone, and the 
model is integrable in terms of elementary functions. The other cases are represented 
generically by a 7^ 6. We will discuss two interesting particular cases, the Landau- 
Lifshitz model with 6 = 0, and the Fedorov model with a — —b. 



3.3.1 Exact integrable one-parameter model with a = b 

In the framework of this one-parameter model equations (46) and (47) reduce to 
two D'Alembert type equations 



dud,[M2 + Q{u) 9„M3] = 0, 
dud,[Ms + Q{u) d,M2] = 0, 



where 



Q{u) 



p* po 



V2L 



V2 



u3u ' 



(48) 
(49) 

(50) 



Solving equations (48) and (49) we obtain 

M2{u, v) + Q{u) d,M^{u, v) = $2(«) + '^2{v) , (51) 

Ms{u, v) + Q{u) d.M^iu, v) = ^s{u) + ^3{v) , (52) 

where $a(ii) and ^a('f) ^-re arbitrary functions of their arguments. The standard 
procedure of the integration of such a system is the following. First, one substitutes 
the M2{u, v) function given in equation (51) into equation (52), and obtains a second 
order differential equation for M^{u,v), 

Ms{u, v) - Q\u) dlM^{u, v) = ^^{u) + *3(^^) - Q{u) d,^2{v) . (53) 

Since the coefficients in (53) depend on u only, one can use the Fourier transforms 
of Ma{u,v) and other functions defined below, in order to find the solutions of the 
respective equation. The relation between a Fourier transform quantity T{u, k) and 
the function T{u, v) is given by 

/OO (jh 
_^-r{u,k)e-^''r (54) 

Applying a Fourier transformation to the quantity M^^u, v), one obtains from equa- 
tion (53), for the Fourier transform M.3{u, k), the formula 

where 

T{u, k) = / dve'^'^ [$3(«) + *3(^) - Q{.u) d,^2{v)] . (56) 

J — OO 

This is valid for arbitrary initial functions ^a{u) and ^a(i')- Let us consider now a 
special set of initial functions, which yield an exact solution in terms of elementary 
functions. We choose ^^(t') as 

*2(v) = *2(0)cos(fcv + (^o), *3(^^) = *3(0)sin(/cv + (/3o) ■ (57) 

We also search for Ma{u, v) in the form 

M^{u, v) = + Xa{u) -^^{v) . (58) 

The functions X^iu) have to be found from the equations (51) and (52). These 
equations yield 



Q{U)X^(U) «2(t.) *2(0) • 

and 



1-X3(.) _%(V)_ ^^^^ 



Q{u)X2{u) -^siv) *3(0) 
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where the right hand side in equations (59) and (60) were put equal to a constant 
of separation of variables, as usual in partial differential equations. These constants 
were chosen in accord the initial conditions given in equation (57). The solutions to 
Xa{u) can be obtained from (59) and (60) and have the form 



^2(0) - kQ{u)^3{0) 



*2(0)(1 + A;2g2(^,)) ' -V ^ ■^s{0){l + k^Q^{u)) 
Thus, from (45), (58), and (61), we obtain for the potentials 



(61) 



^2 = e^ 



A, 



■ ^2(0)-A:g(^)^3(0) ^ ; 

Mu) + i + fe2Q2(,) {kv + ^0) 

^ , , ^3(0)+A;g(ii)*2(0) . ; 
$3(m) + — : , ,^ ' / , sm {kv + ifo) 



(62) 



(63) 



1 + A;2g2(^i) 

The electromagnetic field represented by the formulae (62) and (63) consists of a 
mixture of electromagnetic waves traveling in the same direction and in opposite 
direction to the gravitational pp-wavc direction. 

For ^a{'v) = 0, the electromagnetic wave is a function of the retarded time u 
only, and therefore travels in the direction of the pp-wave. In this case there is no 
optical activity (see, (62) and (63)), since the electromagnetic wave does not couple 
with the curvature (see, also, [18, 19, 20]). 

Thus, let us study the case of electromagnetic waves traveling in an opposite 
direction to the gravitational wave. First, we assume 



^2{u) = ^3{u) = 0. 



(64) 



Then, eliminating the variable v in equations (62) and (63) one obtains the equation 
for the ellipse of polarization 



(65) 



where the quantities 6^X2 (w) and e'^X^iu) play the role of semi-axes for this ellipse. 
The instantaneous position of the polarization vector can be characterized by the 
angle (p{u, k), which satisfy the relationship 



tsinip{u, k) 



X^{u) _ ^3(0) + kQ{u)^<2{Q) 
X2{u) *2(0) - A;g(u)*3(0) 



(66) 



We see, that the polarization vector is rotating. The frequency of rotation, v{u, k), 
is equal to 



i/fu, k) = -^iplu, k) 
as 



d , J . du 
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Thus, 

^(^u,k) = ^^^^ . (68) 

The initial values ^2(0) and ^3(0) are arbitrary. Defining ojcm as the initial frequency 
of the electromagnetic wave, one has k = cUem/c From (50), Q{u) is proportional 
to the Riemann tensor component. Thus, from equation (68), we find that there is 
rotation of the polarization vector. This means there is a net Faraday effect, induced 
by curvature in this instance. 

3.3.2 Solutions for a b 

For a 7^ 6 there are many possible models. We select two. 
(i) The Landau-Lifshitz model: 6=0 

In contrast to the case a — b, the models with a 7^ 6 do not admit a representation of 
the solution in terms of elementary functions. However, for a small enough ratio of 

<^Gw/^cmj where a;Gw is the gravitational wave frequency, one can apply the Fourier 
transformation (54) to the equations (46) and (47). This is a case that occurs in 
astrophysical situations of interest (see, e.g, [25, 26])). Putting 6 = in equations 
(46) and (47) and taking into account the leading terms in k — uJem/c, one obtains 
the following equations for the Fourier transforms A42{u, k) and M.2,{u, k) 

duM2{u, k) = -^Ksu Ms{u, k) , duM2{u, k) = -^Ksu M2{u, k) . (69) 

The solution of (69) is 

M2{u, k) = Al2(0, k) cosh e(it, k) + A^3(0, k) sinh e(w, k) , 

Mi{u, k) = Al3(0, k) cosh Q{u, k) + Al2(0, k) sinh Q{u, k) , (70) 

where 

7 2 

<d{u,k)^^j^duRl,M- (71) 

One sees from (70) and (71) that the plane of polarization of each Fourier component, 
with fixed /c 7^ 0, is rotating, the rotation being hyperbolic. A particular case, 
worth of comment, is when the electromagnetic wave is polarized at m = in the 
x^— direction, and the Fourier amplitude 7W3(0, k) is set to zero, ^^3(0, k) = 0, also 
aiu — 0. Then, from equation (70), one finds that the A43{u, k) component becomes 
nonvanishing at it > 0, since we have assumed A^2(0, k) ^ 0. In this particular case, 
the angle of rotation (/?(■«, A;), defined by tan</?(u, k) = ^^{u, k)/Ai2{u, k), is equal 
to ip{u, k) = arctan (tanh0(M, k)). The corresponding frequency of the rotation of 
the polarization plane, ^{u, k) has the form 

k) = ^ ' / — - . 72 
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(ii) The Fedorov model: a — —b 

In classical electrodynamics the Fedorov model possesses the same duality invari- 
ance, as the Maxwell equations (1). Applying the transformation (see, e.g., [23]) 

E'cos^ + H'sin^, H = -^'sin^ + ^'cos^ , 
D^D'cos^ + B'smC, S = -5'sin^ + S'cos^ (73) 

with arbitrary constant ^ to the Maxwell equations (1), one obtains the same equa- 
tions for the quantities with primes. Analogously, applying the transformations (73) 
to the equations (7) one sees that these equations are duality invariant when a — —b 
only, i.e. in the Fedorov model. 

Let us consider now the curvature generalization of the Fedorov model and solve 
the equations (46) and (47) in the case a = —b. In a leading order approximation in 
k — cuem/c (see part (i) of this subsubsection) one can see that the solutions coincide 
with the solutions describing the Landau-Lifshitz model, but instead of a one has 
to insert 2a. In other words, this model also describes a Faraday rotation of the 
polarization plane with double rotation frequency as comparison to the model 6 = 
(see equation (72)). 



3.3.3 Orders of magnitude for the Faraday rotation 

It is interesting to compare the frequencies of rotation of the polarization vector in 
the model a = 6 on one hand, and the models b — 0, a — —b on the other hand, when 
ai?^„3„ is small. In this case, formula (68) for the model a — b yields a frequency 

1^.=^ - IkaiR'^M)' ■ (74) 

For the Landau and Lifshitz model formula (72) gives 

X 

i/b=o = -k'^aR^^-iM > (75) 
and for the Fedorov model one finds 

i^a=-b = 2 i/b=0 ■ (76) 

One sees that the ratios z/b=o/t'a=b and z/a=-b/t'a=b are of the same order of magnitude 
as the ratio of c<Jcm/^GW- In other words, for cJcm/^GW >> 1) the rate of rotation of 
the polarization vector induced by curvature, predicted in the Landau-Lifshitz and 
Fedorov models, is much faster than the rotation rate predicted in the model with 
a — b. 

Let us estimate numerically the typical magnitude of the effect of the rotation 

of the polarization vector induced by the curvature on the example of the Fedorov 
model. One obtains from (71), (33) and (76) that in the weak gravity field approx- 
imation 

e{u*,k) = ^\P\u*)-(3'm . (77) 
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Consider now the simple model in which the electromagnetic signal emitted by some 
astrophysical source passes through the environment of the pulsar PSR J0534+2200 
in the CRAB nebula and arrives at the Earth. The parameters of this pulsar are 
well-known (e.g. [27]): the period of rotation is T^t = 33.4 ms, the frequency of the 
gravitational radiation is uqw — 60 Hz, the distance is equal to d — 6.2 x 10^^ cm, 
and the optimistic estimation of the amplitude of the gravitational wave on the 
Earth surface is h = 9.5 x 10~^^. 

One can consider the moment of the retarded time u = u* as the moment of the 
arrival of the signal on Earth. As for u = 0, it can be considered as the moment 
when the electromagnetic signal passes in the vicinity of the pulsar. Taking into 
account that the amplitude of the gravitational wave decreases reciprocally to the 
distance from the source, and the gravitational wave field is modeled to be periodic 
with the frequency o^gw — Stt^gW) one can represent the leading term in (77) in the 
following form 

|e(K*)H^^^^^V^, (78) 

where Pq is the amplitude of the gravitational wave at the minimal distance from 
the pulsar (it is of the order of the neutron star radius Rq — 10^ cm), and kc — cUem- 
Putting 

/?o = /i -|- ~ 5.9 X 10"^ , (79) 
Rq 

we obtain the following estimate for the angle of rotation of the polarization vector 
induced by curvature 

|eM|.«.ad(^)(^)(,^)^(|^)^ (80) 

In order to estimate the parameter a one can mention that for a crystal of quartz 
the standard gyration coefficient 7^^^ (see (2)) for the frequency Vem = 5 x 10^^ Hz 
is equal to 7^^^ — 1.2 x 10~^°cm (see, e.g. [28]). This means that the ratio ^quartz = 

123 

^ — , where Aem is the wavelength of the electromagnetic wave, is of the order of 
10~^ . Analogously, we can introduce a new dimensionless parameter ^qw using, for 
example, the following relation 

a = Cgw Aem Xlw - Cgw X 10^^ cm^ , (81) 

where Agw = 7^ ~ 5 x 10^ cm is the wavelength of gravitational radiation. Thus, 
the estimation for the angle rotation yields 

|e(^/*)| ~41radx (10^ x^Gw), (82) 
where {gw << Cquartz < 10~^ has to be estimated experimentally. 
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4 Conclusions 



We have established a model with two gyration parameters a and b, for optical 
activity induced by curvature. Within this two-parameter model we have studied 
three models, namely, a — b, b — 0, and a — —b. It was shown that the a — 
b model is exactly integrable, while for the models 6 = (the Landau-Lifshitz 
model) and a = —b (the Fedorov model) the solutions are presented in the short 
wavelength approximation. The solutions show that the optical activity induced by 
curvature leads to Faraday rotation, i.e., the rotation of the plane of polarization 
of the electromagnetic wave traveling in the gravitational wave background. The 
frequency of the rotation is deternd. The frequency of the rotation is determined by 
the value of the Riemann tensor and by the phenomenological gyration parameters 
a and b. 
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